Multiscale stiffness of human emphysematous precision cut lung slices.

Emphysema is a debilitating disease that remodels the lung leading to reduced tissue stiffness. Thus, understanding emphysema progression requires assessing lung stiffness at both the tissue and alveolar scales. Here, we introduce an approach to determine multiscale tissue stiffness and apply it to precision-cut lung slices (PCLS). First, we established a framework for measuring stiffness of thin, disk-like samples. We then designed a device to verify this concept and validated its measuring capabilities using known samples. Next, we compared healthy and emphysematous human PCLS and found that the latter was 50% softer. Through computational network modeling, we discovered that this reduced macroscopic tissue stiffness was due to both microscopic septal wall remodeling and structural deterioration. Lastly, through protein expression profiling, we identified a wide spectrum of enzymes that can drive septal wall remodeling, which, together with mechanical forces, lead to rupture and structural deterioration of the emphysematous lung parenchyma.

Predicting alveolar ventilation heterogeneity in pulmonary fibrosis using a non-uniform polyhedral spring network model.

Pulmonary Fibrosis (PF) is a deadly disease that has limited treatment options and is caused by excessive deposition and cross-linking of collagen leading to stiffening of the lung parenchyma. The link between lung structure and function in PF remains poorly understood, although its spatially heterogeneous nature has important implications for alveolar ventilation. Computational models of lung parenchyma utilize uniform arrays of space-filling shapes to represent individual alveoli, but have inherent anisotropy, whereas actual lung tissue is isotropic on average. We developed a novel Voronoi-based 3D spring network model of the lung parenchyma, the Amorphous Network, that exhibits more 2D and 3D similarity to lung geometry than regular polyhedral networks. In contrast to regular networks that show anisotropic force transmission, the structural randomness in the Amorphous Network dissipates this anisotropy with important implications for mechanotransduction. We then added agents to the network that were allowed to carry out a random walk to mimic the migratory behavior of fibroblasts. To model progressive fibrosis, agents were moved around the network and increased the stiffness of springs along their path. Agents migrated at various path lengths until a certain percentage of the network was stiffened. Alveolar ventilation heterogeneity increased with both percent of the network stiffened, and walk length of the agents, until the percolation threshold was reached. The bulk modulus of the network also increased with both percent of network stiffened and path length. This model thus represents a step forward in the creation of physiologically accurate computational models of lung tissue disease.

A Personalized Spring Network Representation of Emphysematous Lungs From CT Images.

Emphysema is a progressive disease characterized by irreversible tissue destruction and airspace enlargement, which manifest as low attenuation area (LAA) on CT images. Previous studies have shown that inflammation, protease imbalance, extracellular matrix remodeling and mechanical forces collectively influence the progression of emphysema. Elastic spring network models incorporating force-based mechanical failure have been applied to investigate the pathogenesis and progression of emphysema. However, these models were general without considering the patient-specific information on lung structure available in CT images. The aim of this work was to develop a novel approach that provides an optimal spring network representation of emphysematous lungs based on the apparent density in CT images, allowing the construction of personalized networks. The proposed method takes into account the size and curvature of LAA clusters on the CT images that correspond to a pre-stressed condition of the lung as opposed to a naïve method that excludes the effects of pre-stress. The main findings of this study are that networks constructed by the new method 1) better preserve LAA cluster sizes and their distribution than the naïve method; and 2) predict different course of emphysema progression compared to the naïve method. We conclude that our new method has the potential to predict patient-specific emphysema progression which needs verification using clinical data.

Inflation instability in the lung: an analytical model of a thick-walled alveolus with wavy fibres under large deformations.

Inflation of hollow elastic structures can become unstable and exhibit a runaway phenomenon if the tension in their walls does not rise rapidly enough with increasing volume. Biological systems avoid such inflation instability for reasons that remain poorly understood. This is best exemplified by the lung, which inflates over its functional volume range without instability. The goal of this study was to determine how the constituents of lung parenchyma determine tissue stresses that protect alveoli from instability-related overdistension during inflation. We present an analytical model of a thick-walled alveolus composed of wavy elastic fibres, and investigate its pressure-volume behaviour under large deformations. Using second-harmonic generation imaging, we found that collagen waviness follows a beta distribution. Using this distribution to fit human pressure-volume curves, we estimated collagen and elastin effective stiffnesses to be 1247 kPa and 18.3 kPa, respectively. Furthermore, we demonstrate that linearly elastic but wavy collagen fibres are sufficient to achieve inflation stability within the physiological pressure range if the alveolar thickness-to-radius ratio is greater than 0.05. Our model thus identifies the constraints on alveolar geometry and collagen waviness required for inflation stability and provides a multiscale link between alveolar pressure and stresses on fibres in healthy and diseased lungs.

Tracking respiratory mechanics around natural breathing rates via variable ventilation.

Measuring respiratory resistance and elastance as a function of time, tidal volume, respiratory rate, and positive end-expiratory pressure can guide mechanical ventilation. However, current measurement techniques are limited since they are assessed intermittently at non-physiological frequencies or involve specialized equipment. To this end, we introduce ZVV, a practical approach to continuously track resistance and elastance during Variable Ventilation (VV), in which frequency and tidal volume vary from breath-to-breath. ZVV segments airway pressure and flow recordings into individual breaths, calculates resistance and elastance for each breath, bins them according to frequency or tidal volume and plots the results against bin means. ZVV's feasibility was assessed clinically in five human patients with acute lung injury, experimentally in five mice ventilated before and after lavage injury, and computationally using a viscoelastic respiratory model. ZVV provided continuous measurements in both settings, while the computational study revealed <2% estimation errors. Our findings support ZVV as a feasible technique to assess respiratory mechanics under physiological conditions. Additionally, in humans, ZVV detected a decrease in resistance and elastance with time by 12.8% and 6.2%, respectively, suggesting that VV can improve lung recruitment in some patients and can therefore potentially serve both as a dual diagnostic and therapeutic tool.

Comparison of a micro-electro-mechanical system airflow sensor with the pneumotach in the forced oscillation technique.

PURPOSE: This study supports the use of thin-film micro-electro-mechanical system (MEMS) airflow sensors in the forced oscillation technique. MATERIALS AND METHODS: The study employed static testing using air flow standards and computer-controlled sound attenuations at 8 Hz. Human feasibility studies were conducted with a testing apparatus consisting of a pneumotach and thin-film MEMS air flow sensors in series. Short-time Fourier transform spectra were obtained using SIGVIEW software. RESULTS: Three tests were performed, and excellent correlations were observed between the probes. The thin-film MEMS probe showed superior sensitivity to higher frequencies up to 200 Hz. CONCLUSION: The results suggest that lower-cost thin-film MEMS can be used for forced oscillation technique applications (including home care devices) that will benefit patients suffering from pulmonary diseases such as asthma, COPD, and cystic fibrosis.

The role of heterogeneity in asthma: a structure-to-function perspective.

A number of methods have evolved through the years in probing the dysfunction that impacts mechanics and ventilation in asthma. What has been consistently found is the notion of heterogeneity that is not only captured in the frequency dependence of lung mechanics measurements but also rendered on imaging as patchy diffuse areas of ventilation defects. The degree of heterogeneity has been linked to airway hyperresponsiveness, a hallmark feature of asthma. How these heterogeneous constriction patterns lead to functional impairment in asthma have only been recently explored using computational airway tree models. By synthesizing measurements of lung mechanics and advances in imaging, computational airway tree models serve as a powerful engine to accelerate our understanding of the physiologic changes that occur in asthma. This review will be focused on the current state of investigational work on the role of heterogeneity in asthma, specifically exploring the structural and functional relationships.

Hyperresponsiveness: Relating the Intact Airway to the Whole Lung.

We relate changes of the airway wall to the response of the intact airway and the whole lung. We address how mechanical conditions and specific structural changes for an airway contribute to hyperresponsiveness resistant to deep inspiration. This review conveys that the origins of hyperresponsiveness do not devolve into an abnormality at single structural level but require examination of the complex interplay of all the parts.

Spatial distribution of airway wall displacements during breathing and bronchoconstriction measured by ultrasound elastography using finite element image registration.

With every breath, the airways within the lungs are strained. This periodic stretching is thought to play an important role in determining airway caliber in health and disease. Particularly, deep breaths can mitigate excessive airway narrowing in healthy subjects, but this beneficial effect is absent in asthmatics, perhaps due to an inability to stretch the airway smooth muscle (ASM) embedded within an airway wall. The heterogeneous composition throughout an airway wall likely modulates the strain felt by the ASM but the magnitude of ASM strain is difficult to measure directly. In this study, we optimized a finite element image registration method to measure the spatial distribution of displacements and strains throughout an airway wall during pressure inflation within the physiological breathing range before and after induced narrowing with acetylcholine (ACh). The method was shown to be repeatable, and displacements estimated from different image sequences of the same deformation agreed to within 5.3µm (0.77%). We found the magnitude and spatial distribution of displacements were radially and longitudinally heterogeneous. The region in the middle layer of the airway experienced the largest radial strain due to a transmural pressure (Ptm) increase simulating tidal breathing and a deep inspiration (DI), while the region containing the ASM (i.e., closest to the lumen) strained least. During induced narrowing with ACh, we observed temporal longitudinal heterogeneity of the airway wall. After constriction, the displacements and strain are much smaller than the relaxed airway and the pattern of strains changed, suggesting the airway stiffened heterogeneously.

Linking Ventilation Heterogeneity Quantified via Hyperpolarized 3He MRI to Dynamic Lung Mechanics and Airway Hyperresponsiveness.

Advancements in hyperpolarized helium-3 MRI (HP 3He-MRI) have introduced the ability to render and quantify ventilation patterns throughout the anatomic regions of the lung. The goal of this study was to establish how ventilation heterogeneity relates to the dynamic changes in mechanical lung function and airway hyperresponsiveness in asthmatic subjects. In four healthy and nine mild-to-moderate asthmatic subjects, we measured dynamic lung resistance and lung elastance from 0.1 to 8 Hz via a broadband ventilation waveform technique. We quantified ventilation heterogeneity using a recently developed coefficient of variation method from HP 3He-MRI imaging. Dynamic lung mechanics and imaging were performed at baseline, post-challenge, and after a series of five deep inspirations. AHR was measured via the concentration of agonist that elicits a 20% decrease in the subject's forced expiratory volume in one second compared to baseline (PC20) dose. The ventilation coefficient of variation was correlated to low-frequency lung resistance (R = 0.647, P < 0.0001), the difference between high and low frequency lung resistance (R = 0.668, P < 0.0001), and low-frequency lung elastance (R = 0.547, P = 0.0003). In asthmatic subjects with PC20 values <25 mg/mL, the coefficient of variation at baseline exhibited a strong negative trend (R = -0.798, P = 0.02) to PC20 dose. Our findings were consistent with the notion of peripheral rather than central involvement of ventilation heterogeneity. Also, the degree of AHR appears to be dependent on the degree to which baseline airway constriction creates baseline ventilation heterogeneity. HP 3He-MRI imaging may be a powerful predictor of the degree of AHR and in tracking the efficacy of therapy.

Can breathing-like pressure oscillations reverse or prevent narrowing of small intact airways?

Periodic length fluctuations of airway smooth muscle during breathing are thought to modulate airway responsiveness in vivo. Recent animal and human intact airway studies have shown that pressure fluctuations simulating breathing can only marginally reverse airway narrowing and are ineffective at protecting against future narrowing. However, these previous studies were performed on relatively large (>5 mm diameter) airways, which are inherently stiffer than smaller airways for which a preponderance of airway constriction in asthma likely occurs. The goal of this study was to determine the effectiveness of breathing-like transmural pressure oscillations to reverse induced narrowing and/or protect against future narrowing of smaller, more compliant intact airways. We constricted smaller (luminal diameter = 2.92 ± 0.29 mm) intact airway segments twice with ACh (10(-6) M), once while applying tidal-like pressure oscillations (5-15 cmH2O) before, during, and after inducing constriction (Pre + Post) and again while only imposing the tidal-like pressure oscillation after induced constriction (Post Only). Smaller airways were 128% more compliant than previously studied larger airways. This increased compliance translated into 196% more strain and 76% greater recovery (41 vs. 23%) because of tidal-like pressure oscillations. Larger pressure oscillations (5-25 cmH2O) caused more recovery (77.5 ± 16.5%). However, pressure oscillations applied before and during constriction resulted in the same steady-state diameter as when pressure oscillations were only applied after constriction. These data show that reduced straining of the airways before a challenge likely does not contribute to the emergence of airway hyperreactivity observed in asthma but may serve to sustain a given level of constriction.

Correlated variability in the breathing pattern and end-expiratory lung volumes in conscious humans.

In order to characterize the variability and correlation properties of spontaneous breathing in humans, the breathing pattern of 16 seated healthy subjects was studied during 40 min of quiet breathing using opto-electronic plethysmography, a contactless technology that measures total and compartmental chest wall volumes without interfering with the subjects breathing. From these signals, tidal volume (VT), respiratory time (TTOT) and the other breathing pattern parameters were computed breath-by-breath together with the end-expiratory total and compartmental (pulmonary rib cage and abdomen) chest wall volume changes. The correlation properties of these variables were quantified by detrended fluctuation analysis, computing the scaling exponenta. VT, TTOT and the other breathing pattern variables showed α values between 0.60 (for minute ventilation) to 0.71 (for respiratory rate), all significantly lower than the ones obtained for end-expiratory volumes, that ranged between 1.05 (for rib cage) and 1.13 (for abdomen) with no significant differences between compartments. The much stronger long-range correlations of the end expiratory volumes were interpreted by a neuromechanical network model consisting of five neuron groups in the brain respiratory center coupled with the mechanical properties of the respiratory system modeled as a simple Kelvin body. The model-based α for VT is 0.57, similar to the experimental data. While the α for TTOT was slightly lower than the experimental values, the model correctly predicted α for end-expiratory lung volumes (1.045). In conclusion, we propose that the correlations in the timing and amplitude of the physiological variables originate from the brain with the exception of end-expiratory lung volume, which shows the strongest correlations largely due to the contribution of the viscoelastic properties of the tissues. This cycle-by-cycle variability may have a significant impact on the functioning of adherent cells in the respiratory system.

Tidal stretches differently regulate the contractile and cytoskeletal elements in intact airways.

Recent reports suggest that tidal stretches do not cause significant and sustainable dilation of constricted intact airways ex vivo. To better understand the underlying mechanisms, we aimed to map the physiological stretch-induced molecular changes related to cytoskeletal (CSK) structure and contractile force generation through integrin receptors. Using ultrasound, we measured airway constriction in isolated intact airways during 90 minutes of static transmural pressure (Ptm) of 7.5 cmH2O or dynamic variations between Ptm of 5 and 10 cmH20 mimicking breathing. Integrin and focal adhesion kinase activity increased during Ptm oscillations which was further amplified during constriction. While Ptm oscillations reduced ß-actin and F-actin formation implying lower CSK stiffness, it did not affect tubulin. However, constriction was amplified when the microtubule structure was disassembled. Without constriction, α-smooth muscle actin (ASMA) level was higher and smooth muscle myosin heavy chain 2 was lower during Ptm oscillations. Alternatively, during constriction, overall molecular motor activity was enhanced by Ptm oscillations, but ASMA level became lower. Thus, ASMA and motor protein levels change in opposite directions due to stretch and contraction maintaining similar airway constriction levels during static and dynamic Ptm. We conclude that physiological Ptm variations affect cellular processes in intact airways with constriction determined by the balance among contractile and CSK molecules and structure.

A computational model of the response of adherent cells to stretch and changes in substrate stiffness.

Cells in the body exist in a dynamic mechanical environment where they are subject to mechanical stretch as well as changes in composition and stiffness of the underlying extracellular matrix (ECM). However, the underlying mechanisms by which cells sense and adapt to their dynamic mechanical environment, in particular to stretch, are not well understood. In this study, we hypothesized that emergent phenomena at the level of the actin network arising from active structural rearrangements driven by nonmuscle myosin II molecular motors play a major role in the cellular response to both stretch and changes in ECM stiffness. To test this hypothesis, we introduce a simple network model of actin-myosin interactions that links active self-organization of the actin network to the stiffness of the network and the traction forces generated by the network. We demonstrate that such a network replicates not only the effect of changes in substrate stiffness on cellular traction and stiffness and the dependence of rate of force development by a cell on the stiffness of its substrate, but also explains the physical response of adherent cells to transient and cyclic stretch. Our results provide strong indication that network phenomena governed by the active reorganization of the actin-myosin structure plays an important role in cellular mechanosensing and response to both changes in ECM stiffness and externally applied mechanical stretch.

Can tidal breathing with deep inspirations of intact airways create sustained bronchoprotection or bronchodilation?

Fluctuating forces imposed on the airway smooth muscle due to breathing are believed to regulate hyperresponsiveness in vivo. However, recent animal and human isolated airway studies have shown that typical breathing-sized transmural pressure (Ptm) oscillations around a fixed mean are ineffective at mitigating airway constriction. To help understand this discrepancy, we hypothesized that Ptm oscillations capable of producing the same degree of bronchodilation as observed in airway smooth muscle strip studies requires imposition of strains larger than those expected to occur in vivo. First, we applied increasingly larger amplitude Ptm oscillations to a statically constricted airway from a Ptm simulating normal functional residual capacity of 5 cmH2O. Tidal-like oscillations (5-10 cmH2O) imposed 4.9 ± 2.0% strain and resulted in 11.6 ± 4.8% recovery, while Ptm oscillations simulating a deep inspiration at every breath (5-30 cmH2O) achieved 62.9 ± 12.1% recovery. These same Ptm oscillations were then applied starting from a Ptm = 1 cmH2O, resulting in approximately double the strain for each oscillation amplitude. When extreme strains were imposed, we observed full recovery. On combining the two data sets, we found a linear relationship between strain and resultant recovery. Finally, we compared the impact of Ptm oscillations before and after constriction to Ptm oscillations applied only after constriction and found that both loading conditions had a similar effect on narrowing. We conclude that, while sufficiently large strains applied to the airway wall are capable of producing substantial bronchodilation, the Ptm oscillations necessary to achieve those strains are not expected to occur in vivo.

A mechanical design principle for tissue structure and function in the airway tree.

With every breath, the dynamically changing mechanical pressures must work in unison with the cells and soft tissue structures of the lung to permit air to efficiently traverse the airway tree and undergo gas exchange in the alveoli. The influence of mechanics on cell and tissue function is becoming apparent, raising the question: how does the airway tree co-exist within its mechanical environment to maintain normal cell function throughout its branching structure of diminishing dimensions? We introduce a new mechanical design principle for the conducting airway tree in which mechanotransduction at the level of cells is driven to orchestrate airway wall structural changes that can best maintain a preferred mechanical microenvironment. To support this principle, we report in vitro radius-transmural pressure relations for a range of airway radii obtained from healthy bovine lungs and model the data using a strain energy function together with a thick-walled cylinder description. From this framework, we estimate circumferential stresses and incremental Young's moduli throughout the airway tree. Our results indicate that the conducting airways consistently operate within a preferred mechanical homeostatic state, termed mechanical homeostasis, that is characterized by a narrow range of circumferential stresses and Young's moduli. This mechanical homeostatic state is maintained for all airways throughout the tree via airway wall dimensional and mechanical relationships. As a consequence, cells within the airway walls throughout the airway tree experience similar oscillatory strains during breathing that are much smaller than previously thought. Finally, we discuss the potential implications of how the maintenance of mechanical homeostasis, while facilitating healthy tissue-level alterations necessary for maturation, may lead to airway wall structural changes capable of chronic asthma.